Optical imaging beyond the diffraction limit via dark states 
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We study the possibility of creating spatial patterns having subwavelength size by using the so- 
called dark states formed by the interaction between atoms and optical fields. These optical fields 
have a specified spatial distribution. Our experiments in Rb vapor display spatial patterns that 
are smaller than the length determined by the diffraction limit of the optical system used in the 
experiment. This approach may have applications to interference lithography and might be used in 
coherent Raman spectroscopy to create patterns with subwavelength spatial resolution. 
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The ability to create small images is important for ma- 
terial processing technology and for improving the reso- 
lution of microscopy for bio- medical applications [l| . Re- 
cently, several methods have been presented that are able 
to overcome the diffraction limit of the imaging system. 
Quantum microscopy is based on using a nonclassical op- 
tical field approach [j, 0] • Microscopy with classical fields 
can be enhanced by the nonlinear optical response of the 
medium [1] . Classical field amplitude and phase arrange- 
ments can be used to locate the position of an atom with 
subwavelength precision in an atomic beam 0, @, 0, @] 
and in a cavity [9j]. 

Here, we suggest a new approach that is based on co- 
herent population trapping [TO, EH, EH, E, 03 • Optical 
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FIG. 1: Qualitative description of the idea, (a) Distribution 
of the drive (1) and the probe (2) fields vs a transverse spatial 
coordinate at the entrance to the cell, (b) Dependance of the 
absorption coefficient given by Eq.j4} vs position. Plots (c) 
and (d) show the distribution of the probe beam after prop- 
agating through the cell. Case (c) is for a strong drive field 
and relatively low optical density. Case (d) is for a relatively 
weak drive field and large optical density. 



fields applied to a three-level quantum system excite the 
so-called dark state, which is decoupled from the fields. 
Similar approaches using coherent population trapping 
have also been developed by several groups (for example, 

see [IE EE 113) • 

As a qualitative introduction, assume that the drive 
field Rabi frequency fid has the particular spatial distri- 
bution sketched in Fig. [2(a) by the solid line (1). The 
weak probe field Rabi frequency fl p (fl p <C fid) has a 
diffraction limited distribution (shown by the dashed line 
(2) in Fig. [T](a)). The probe and drive fields are applied 
to the atom (see the inset in Fig. [2]). At all positions 
of nonzero drive field, the dark state, which is given [l2[ 



by \D) = (fl p \c) - fl d \b))/^fll + fl% is practically \b). 

When the drive field is zero, the dark state is |c), and 
the atoms at these positions are coupled to the fields and 
some atoms are in the upper state \a). The size of a spot 
where the atoms are excited depends on the relaxation 
rate 7 C & between levels \b) and \c). For >y cb = 0, the size 
of spot is zero, smaller than the optical wavelength. 

The Hamiltonian of three-level atom interacting with 
optical fields (see the inset in Fig. [2]) is given by 



H = Ml d \a){b\ + hft p \a)(c\ + adj., 



(1) 



where f2 p>t / = p v ^E v ^jh are the Rabi frequencies of the 
drive Ed and the probe E p fields, respectively. Then, 
the atomic response is given by the set of density matrix 
equations [12| 
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where F describes the relaxation processes. The propa- 
gation of the probe field fl p through the cell is governed 
by Maxwell's Equations and, for propagation in the z- 
direction, can be written in terms of the probe field Rabi 
frequency as 



dftp 

dz 
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The first term accounts for the dispersion and absorption 
of the resonant three-level medium, and the second term 
describes the focusing and/or diffraction of the probe 
beam. The density matrix element p a b is related to the 
probe field absorption which in turn depends on the de- 
tuning and the drive field. This is characterized by an 
absorption coefficient: 
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where T C 5 = 7 C ^ + iuo and T a 5 = 7 + ico; uo = oo^ — 
v is the detuning from the atomic frequency uj a b; 7 is 
the relaxation rate at the optical transition; and 77 = 
3A 2 A^7 r /87r; j r is the spontaneous emission rate. We now 
assume that the drive field has a distribution of intensity 
near its extrema given by 



\n d (z,x)\ 2 = \n f 
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where = Qd( z , x o), a typical absorption profile vs x 
is shown in Fig. []Jb). Neglecting the diffraction term in 
Eq.Q, we can write an approximate solution for Eq. ([3|) 

as 



£l p (z, x) = Q p (z = 0, x) ex.p(—Kz). 



(6) 



For relatively low optical density (nz ~ 1), nearly all 
of the probe field propagates through the cell except for 
a small part where the drive field is zero (see Fig. Ufa)). 
Absorption occurs there because the probe beam excites 
the atomic medium. The width of the region of the ex- 
cited medium, in the vicinity of zero drive field, is char- 
acterized by 
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where ft = Qd(z = 0, x = 0). This region is small, but 
its contrast is limited because of the finite absorption of 
the medium at the center of optical line (Fig. [1(c)). 

For higher optical density, this narrow feature becomes 
broadened (compare Fig. [Tfc and d), but two narrow 
peaks are formed during the propagation of the probe 
beam (see Fig. [I^d)). For zero detuning, their width is 
given by 



Ax = L 
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FIG. 2: (Color online) Experimental schematic. A/2: half- 
wave plate; A/4: quarter- wave plate; LI, L2, L3: lenses; 
MZ: Mach-Zehnder interferometer; PZT: piezo-electric ele- 
ment; PBS: polarizing beam splitter, PD: photo diode; CCD: 
CCD camera. Picture A is the spatial intensity distribution 
of the drive field. Picture B is the beam profile of the paral- 
lel probe beam without the lens LI. Picture C is the beam 
profile of the diffraction limited probe beam with the lens 
LI. The inset is the energy diagram of the Rb atom, showing 
representative sublevels. 
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FIG. 3: (Color online) The results of the experiment with 
a parallel probe beam. Picture (a) shows the image of the 
intensity distribution of the drive field in the Rb cell. Picture 
(b) shows the intensity distribution of the transmitted probe 
field. Curves (c) and (d) are the corresponding profiles. The 
widths of the peaks in curves (c) and (d) are 0.4 mm and 0.1 
mm, respectively. 



Note that the propagation occurs in a waveguide 
formed by the coherent medium. The drive field provides 
flexibility for creating patterns with sizes smaller than 
the wavelength of the laser. The distribution of fields is 
governed by electrodynamics and has a diffraction limit, 
while the distribution of molecules in their excited states 



is NOT related to the diffraction limit, but rather deter- 
mined by the relaxation rates T a 5 and T C 5, and thus can 
have spatial sizes smaller than the wavelength. 

In this Letter, we report a proof-of-principle experi- 
ment in Rb vapor to demonstrate our approach. We have 
observed that the distribution of the transmitted probe 
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beam intensity has a double-peak pattern, which is sim- 
ilar to that of the drive beam, but the width of peaks of 
the probe beam is narrower than that of the drive beam. 

The experimental schematic is shown in Fig. O To ob- 
tain diffraction limited focusing the laser beam must have 
a good quality spatial profile. This is obtained by send- 
ing the radiation of an external cavity diode laser through 
a polarization-preserving single mode optical fiber. The 
laser beam is vertically polarized and split into two beams 
(drive and probe). The probe beam carries a small por- 
tion of the laser intensity, and its polarization is rotated 
to be horizontal. 

To create a double-peak spatial distribution for the 
drive field, the drive beam is split into two beams that 
cross at a small angle, using a Mach-Zehnder interferom- 
eter (shown in the dashed square of Fig. [2j). A typical 
two-peak interference pattern of crossing beams is shown 
as Fig. [2]A. 

The probe and drive beams combine on a polarizing 
beam splitter, arranged so that the probe field and the 
interference pattern of the drive field are overlapped in a 
Rb cell. The Rb cell has a length of 4 cm, and is filled 
with 87 Rb. A magnetic shield is used to isolate the cell 
from any environmental magnetic fields, while a solenoid 
provides an adjustable, longitude magnetic field. The 
cell is installed in an oven that heats the cell to reach 
an atomic density of 10 12 cm -3 . The laser is tuned to 
the Di line of 87 Rb at the transition 5 2 S 1/2 (F = 2) -> 
5 2 P 1/2 (F = 1). 

As stated above, the probe and drive beams have the 
orthogonal linear polarizations. A quarter- wave plate 
converts them into left and right circularly polarized 
beams, which couple two Zeeman sublevels of the lower 
level and one sublevel of the excited level of the Rb atoms 
(see the inset of Fig. [2j). 

After passing through the cell, the probe and drive 
beams are converted back to linear polarizations by an- 
other quarter- wave plate and the separated by a polariz- 
ing beam splitter (PBS). The power of transmitted probe 
field is monitored by a photodiode (PD). The spatial in- 
tensity distribution of probe field is recorded by an imag- 
ing system, consisting of the lens L3 and a CCD camera. 

The intensity of the probe beam is low enough that its 
transmission through the cell is almost zero without the 
presence of drive laser. Applying the drive laser makes 
the atomic medium transparent for the probe laser wher- 
ever the EIT condition is satisfied. If the drive laser has a 
certain transverse spatial distribution, then that pattern 
can be projected to the transmission profile of the probe 
laser. 

Two different experiments have been performed. In 
the first experiment, the lenses LI and L2 are not used, 
and the probe beam is a parallel beam with a diameter 
of 1.4 mm. The image of the drive intensity distribution 
in the cell is shown in Fig.[3](a). The probe intensity has 
a Gaussian distribution before entering the cell, and its 
distribution is similar to the drive intensity distribution 
after the cell. As shown in Fig. 09(b), however, the trans- 



mitted probe intensity has a distribution that has sharper 
peaks compared with the pattern of the drive intensity. 
The horizontal cross-sections of the drive and the trans- 
mitted probe distributions are shown in Fig. [31(c) and 
(d) respectively. In the drive intensity profile, the width 
(FWHM) of the peaks is 0.4 mm. The width (FWHM) 
of the peaks in the transmitted probe intensity profile is 
0.1 mm. The spacing between two peaks are the same for 
both the drive and transmitted probe fields. We define 
the finesse as the ratio of the spacing between peaks to 
the width of peaks. The finesse of the transmitted probe 
intensity distribution is a factor of 4 smaller than that of 
the drive intensity distribution. 
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FIG. 4: (Color online) The results of the experiment with 
the diffraction limited probe beam. Picture (a) shows the 
image of the intensity distribution of the drive field in the Rb 
cell. Picture (b) shows the image of the intensity distribution 
of the transmitted probe field. Curves (c) and (d) are the 
corresponding profiles. The widths of the peaks in curves (c) 
and (d) are 165 /xm and 93 /im, respectively. 

In the second experiment, the lenses LI and L2 are 
used. A parallel probe beam (Fig. [2f>) with a diameter 
of 1.4 mm is focused by the lens LI, which has a focal 
length of 750 mm. The beam size at the waist is 0.5 mm, 
which is diffraction limited. To assure experimentally 
that the beam is diffraction limited, we increased the 
beam diameter of the parallel beam by the factor of 2, 
and the beam size at the waist became two times smaller. 
The lens L2 is used to make the drive beam smaller (but 
not diffraction limited) in the Rb cell, where the pattern 
of drive field is spatially overlapped with the waist of the 
probe beam. Classically, there should be no structures at 
the waist of the probe beam because it is diffraction lim- 
ited. Structures can be created in a region smaller than 
the diffraction limit in our experiment, however. The 
experimental result is shown in Fig. 2J The drive field 
still has a double peak intensity distribution (Fig. Ufa)). 
The transmission of the diffraction limited probe beam 
also has a double-peak intensity distribution as shown 
in Fig. 2Kb). Curves (c) and (d) are the beam profiles 
of the drive and transmitted probe beams respectively. 
The width of the peaks in the drive beam is 165 /im, and 
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FIG. 5: Narrowing of the transmitted probe intensity distri- 
bution as function of the probe detuning: (a) experimental 
results (a) and (b) theoretical simulation. The transmitted 
probe profile is shown as well. 

the width of the peaks in the transmitted probe beam is 
93 /im. The finesse of the transmitted probe beam is 1.8 
times greater that that of the drive beam. For the probe 
beam, the structure created within the diffraction limit 
has a size characterized by the width of peaks (93 /im). 
This characteristic size is 5 times smaller than the size of 
the diffraction limited probe beam (500 /im). 

Thus, we have demonstrated that our concept works in 
Rb vapor. The typical image is shown in Fig.[2j One can 
see that the width of the probe image (C) is at least three 
times smaller than the width of the drive image (A). Al- 
though the diffraction limit is "beaten," the experiment 
does not violate any laws of optics. The probe beam is 
diffraction limited, but the atoms are much smaller than 
the size of diffraction limited beam. Moreover, due to 



the strong nonlinearity of the EIT, the characteristic size 
of the pattern in the transmitted probe beam is much 
smaller than that of the drive beam and the diffraction 
limit of the probe beam. 

We have also measured the narrowing effect vs the de- 
tuning of the probe field and have performed simulations 
using the density matrix approach. The results are shown 
in Fig. [5l The calculations reproduce the data satisfac- 
torily. The dependence on detuning has not been consid- 
ered in [IEEEEzlEi|. it is unique for our approach and 
can be understood in the following way. Absorption by 
the atomic medium given by Eq. (|4j) with a drive intensity 
distribution given by Eq.© can be written as 

Then, ratio of the width of the probe intensity distribu- 
tion to the width of the drive intensity distribution is 
given by 

From this we see that the finesse increases with the de- 
tuning. 

In conclusion, we have performed a proof-of-principle 
experiment that our concept works in Rb vapor and have 
experimentally demonstrated the possibility of creating 
structures having widths smaller than those determined 
by the diffraction limits of the optical systems. The re- 
sults obtained here can be viewed as an experimental 
verification of our approach, as well as evidence support- 
ing the theoretical predictions and results obtained by 
others [HI, 

EE E3, El- This technique might be used 
in microscopy by studying the distribution of molecules 
with subwavelength resolution or in lithography by ma- 
nipulating molecules in the excited state. Also, note that 
it may be possible to apply this approach to coherent Ra- 
man scattering (for example, CARS). This may improve 
the spatial resolution of CARS microscopy. 
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